DTP/96/44 
RAL-TR-96-037 
May 1996 



<o '. Parton distributions: a study of the new 

HERA data, a s , the gluon and pp jet production 



A. D. Martin , R. G. Roberts 6 and W. J. Stirling a - C 

a Department of Physics, University of Durham, Durham, DH1 3LE 
b Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX 
i/^ \ c Department of Mathematical Sciences, University of Durham, Durham, DH1 3LE 

m 
\o 
o 
\o 

a: 



X 



Abstract 

New data, especially HERA measurements of the proton structure function at small 
x, allow the opportunity to improve our knowledge of the gluon and quark distribution 
functions. We perform a global analysis which incorporates these new precise data, and 
which extends down to Q 2 = 1.5 GeV 2 . We discuss the sensitivity to the value of a s and 
the improvement in the determination of the gluon. We compare the predictions of the 
single jet inclusive cross section with recent measurements from the Fermilab Tevatron. 



Predictions for hard scattering processes involving hadrons rely on a precise knowledge of 
the parton distribution functions. Traditionally these are determined from global analyses of 
data for a wide range of processes, see for example JT], |2|. As the precision of data improves 
and the kinematic domain of the measurements enlarges we obtain a more rigorous test of 
perturbative QCD and a better knowledge of the parton distributions. This is well illustrated 
by the new measurements |3|, |J of the proton structure function at HERA, which are more 
precise and extend to lower x than hitherto. 

The significant improvement in the data demands a concomitant refinement of the parton 
distributions. In this Letter we describe a new global analysis, based on next-to-leading-order 
DGLAP evolution, of the type described in Ref. [|IJ but with the following special features: 

(i) structure function data down to Q 2 = 1.5 GeV 2 are included; 

(ii) fits are performed for two different values of the QCD coupling, the standard "deep 
inelastic scattering" value (a s (M|) = 0.113) [[J and a larger value (a s (M§) = 0.120) 
suggested, for example, by measurements at LEP 0; 

(iii) particular attention is paid to the small x behaviour of the gluon and sea quark distribu- 
tions; 

(iv) the single jet inclusive cross section is predicted and compared with the new CDF f?J and 
DO || measurements at the Fermilab pp collider. 



We elaborate on these aspects of the analysis in turn below. 

In previous analyses we fitted to structure function data with Q 2 > 5 GeV 2 .[] However 
the "dynamical" partons of the GRV model |J were found to give a good description of data 
down to Q 2 ~ 1 GeV 2 , although the new HERA measurements |3], |J at low Q 2 do show 
some systematic discrepancy between the GRV predictions and the data. In this respect it 
is important to investigate whether this signals a problem for DGLAP evolution at such low 
scales or for the assumptions implicit in the dynamical parton model. We therefore include 
all data with Q 2 > 1.5 GeV 2 . The choice cn s (M§) = 0.113 corresponds to the optimum QCD 
coupling as determined by the scaling violations of the fixed-target deep inelastic data. The 
key data here are the BCDMS _F 2 MP ' d measurements fH| in the interval 0.3 < x < 0.5. There 
is, however, some indication that the HERA structure function data and the Fermilab jet data 



prefer a larger value of the coupling |L1|, [12], [13], [14[] . To illustrate this point we therefore also 
present results of a second global fit in which the coupling is fixed at a s (M|) = 0.120. We call 
the resulting set of partons^ Ri and R2 respectively. 

1 Contamination by higher-twist contributions at large x is avoided by the additional cut W 2 = Q 2 (l — x)/x > 
10 GeV 2 . 

2 These new distributions were presented in preliminary form at the Rome DIS96 Workshop. 
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The advent of the HERA measurements of F 2 has considerably improved our knowledge of 
the gluon and sea quark distributions in the small x regime. We may represent their leading 
small x behaviour by the forms 

xg -> X - X ^ Q2) , xS -> x- Xs(Q2 *> . (1) 

In the perturbative region the evolution of the sea quark distributions is driven by the gluon, 
which is the dominant parton at small x, via the g — > qq splitting. The distributions are 
therefore correlated at small x, and for this reason most previous global analyses set As = X g in 
([]]) at the input scale, usually taken to be Ql = 4 GeV 2 . Even disregarding non-perturbative 
effects, this expectation is only approximate. As we evolve up in Q 2 using the DGLAP evolution 
equations, dg/d\riQ 2 = P gg ® g + . . . and dq/d\nQ 2 = P qg <g> g + . . ., the effective exponents 
\{Q 2 ) increase with Q 2 , and satisfy 

As(Q 2 ) = \ 9 {Q 2 ) ~ e, (2) 

where e is positive and slowly varying with Q 2 , provided that the evolution is sufficiently long. 
Typically we find e ~ 0.05 at Q 2 ~ 200 GeV 2 , see below. That is, DGLAP evolution leads to 
a sea distribution which is slightly less steep in x than the gluon distribution. 

The precision of the new HERA small x structure function data is now such that we can 
begin to study the small x behaviour of the quarks and gluons in detail. We therefore allow 
both As and \ g to be free parameters at the input scale, which we here take to be Ql = 1 GeV 2 . 
To see the extent to which As and X g can be individually determined, we also perform fits with 
As = X g at Qq. Where appropriate, we show properties of these additional parton sets, which 
we label R3 and R4 corresponding to a s (M|) = 0.113 and 0.120 respectively. 

Another interesting development is the constraint that the high statistics jet data from 
Fermilab are starting to impose on the partons and the QCD coupling a s . The slope of the 
inclusive jet transverse energy (E T ) distribution is particularly sensitive to the value of a s [O, 



P^ ]. In the region where the data are most definitive, 50 < Et < 200 GeV, the jets arise from 
gluon-initiated QCD processes. Hence the data can give a tight constraint on a s (fi 2 )g(x, /i 2 ) 
where the gluon is sampled in the region x ~ 2E T /\/s ~ 0.1, but at a high scale /1 2 ~ E\ 
where evolution tends to wash out the differences between the parton sets. We compare the 
predictions based on our new Ri and R2 sets with CDF and preliminary DO jet Et 
distributions to demonstrate the potential discriminating power of the data. 

As in our previous studies, we adopt the following simple analytic form for the starting 
distributions :0 

xfi(x,Q 2 ) = A iX - x *(l - x)^(l + + %x) , (3) 

for i = uy, dy, S, g, but at the lower scale Ql = 1 GeV 2 . Here 5* represents the total sea 
quark distribution. At Ql, only the u, d, s distributions are non-zero, and so 5* = 2(u + d + s). 

3 The distributions arc defined in the MS renormalization and factorization schemes. 
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The flavour structure of the sea is taken to be 



2u = OAS -A 
2d = 0.45 + A 
2s = 0.25, 



with 



xA 



x(d 



u 



0.3/ 



X 



1 + 7a^) 



(4) 



(5) 



The suppression of the strange distribution is motivated by data on neutrino-induced deep 
inelastic dimuon production obtained by the CCFR collaboration [T5J . In our previous studies 



this factor 2 suppression was applied at = 4 GeV , rather than at Qq — 1 GeV as in 
(£|). The consequence is that the strange sea is slightly larger at Q 2 = 4 GeV 2 and, in fact, 
in better agreement with the CCFR measurement. That is, the strange sea lies more centrally 
in the allowed band for s(x,Q 2 = 4 GeV 2 ), determined by the next-to-leading order CCFR 
analysis, than does the MRS (A) curve shown in Fig. 4 of Ref. QXJ . The u, d flavour symmetry 
breaking distribution A(x) is chosen to give good agreement with the NMC measurement of 
the Gottfried sum ||16|| , and the NA51 measurement of the Drell-Yan pp/pn asymmetry fI7 |. 

In previous analyses, we treated the charm quark as massless and evolved its distribution 
from zero at Q z 



m 



where m 2 = 2.7 GeV 2 was chosen to reproduce the EMC charm data 
Now that we have lowered the starting scale to Qo = 1 GeV 2 it is necessary to smooth 
out the onset of charm. We take c(x, Q^) = and evolve up in Q 2 with four flavours but since 
Qq is now much lower, this unmodified density would overestimate the observed F^x^Q 2 ). 
We simply suppress the generated charm density by a smooth factor 
threshold behaviour of the massive quark, i.e. 



Til which mimics the 



where 



c(x, Q z 



Af(z) = 6z 



1 




2z 



' v / 4^TT+ l ' 
vTO ~ V v/4F+T- l 



In 



(6) 



(7) 



with itiq = 3.5 GeV 2 chosen give a reasonable description of the EMC Jl8 and preliminary HI 



20] measurements of the charm structure function F 2 C . We see from Fig. 1 that this approximate 



treatment of charm does indeed give a satisfactory description of F% at present, but as the data 
improve it will become necessary to give a more satisfactory treatment of charm-mass effects. 

Using the above parton parametrizations at Q\ = 1 GeV 2 we perform global next-to-leading 
order DGLAP fits to the deep inelastic and related hard scattering data that were used in the 
MRS(A) analysis flj, but now including data down to Q 2 = 1.5 GeV 2 , supplemented by SLAC 
E665 |E| and updated NMC [2]J structure function measurements and, most importantly, 
the new HI and ZEUS || F 2 data. We present the results of the four different fits Ri - 
R4, which we described above. The optimum values of the starting distribution parameters are 
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listed in Table |l[ Note that three of the four Aj coefficients are determined by the momentum 
and flavour sum rules. The fourth (As) is a fitted parameter. Because of the different values 
of Qq, these parameters cannot be directly compared with the corresponding parameters from 
previous MRS analyses. The \ 2 values found in the four fits for the deep inelastic structure 
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0.113 


0.120 


0.113 
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A ST 4 ( MeV ) 


241 


344 
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-3.85 
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6.55 


6.47 


6.83 


18.1 




At, 


0.60 


0.61 


0.59 


0.66 




Vu 


3.69 


3.54 


3.67 


3.55 






-1.18 


-0.98 


-1.02 


-1.08 


Valence 


In 


6.18 


6.51 


6.36 


5.43 




\l 


0.24 


0.24 


0.25 


0.26 




Vd 


4.43 


4.21 


4.45 


4.18 






5.63 


7.37 


4.82 


9.64 




Id 


25.5 


29.9 


23.5 


26.3 




A s 


0.42 


0.37 


0.92 


0.92 




As 


0.14 


0.15 


0.04 


0.04 




Vs 


9.04 


8.27 


9.38 


8.93 


Sea 


es 


1.11 


1.13 


-1.65 


-2.34 




Is 


15.5 


14.4 


11.8 


12.0 




A A 


.039 


.036 


.040 


.038 




7a 


64.9 


64.9 


64.9 


64.9 



Table 1: The numerical values of the starting distribution parameters of the four sets of partons. 
Note that A g is fixed by the momentum sum rule, and is therefore not a free parameter. 

function subsets of the data are listed in Table 0. We have not included the comparisons of 
each fit with other constraints such as prompt photon production, W asymmetry, di-lepton 
production, and Drell-Yan pp/pn asymmetry, but in every case the quality of the descriptions 
is close to that of MRS (A, A') ||, ||. 

An idea of the quality of the description of the small x measurements of that is obtained 
in the Ri and R2 fits can be seen from the compilation shown in Fig. 2. For comparison we also 
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85 
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8 


SLAC 
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24 


222 


52 


169 


93 


£ T > 200 GeV 


11 


10 


20 


9 


7 
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50<E T < 200 GeV 


14 


103 


72 


50 


126 


E T > 200 GeV 


12 


26 


22 


12 


41 



Table 2: x 2 values for some of the data used in the global fit. The da/dEx jet data are 
predicted, not fitted. 

show the prediction of the GRV set of ("dynamical") partons f9|. In the GRV analysis charm 
is treated as a heavy quark (rather than as a parton), this contribution being included in the 
GRV curves shown. 

Fig. 3 displays the Q 2 dependence of the exponents X g and \s of the x~ x behaviour of 
the gluon and sea quark distributions, see Eq. (|I]). To be precise the curves are obtained by 
assuming that the parton forms @ describe the evolved distributions at any Q 2 and then 
determining the Xi(Q 2 ) by a five-parameter (Ai, Aj, rji, ji) fit to xg(x,Q 2 ) and xS(x, Q 2 ) at 
each Q 2 . The upper plot is obtained from the R 2 set of partons, which have A5 7^ \ g at 
Qq = 1 GeV 2 . It shows that the gluon is valence-like at the input scale. However evolution in 
Q 2 rapidly steepens the small x shape of the gluon and by Q 2 ~ 2 GeV 2 it is already starting 
to drive the sea quark distribution, leading to the DGLAP expectation shown in Eq. (0). It 
is interesting to note that the cross-over point, X g = \s, occurs in the region Q 2 ~ 5 GeV 2 . 
Thus the behaviour is approximately compatible with the MRS(A,A') sets of partons which 
were required to have \ g = A5 at the Qq = 4 GeV 2 starting scale. The data points shown in 
Fig. 3 are obtained by the HI collaboration by fitting their measurements of F^ for x < 0.1 
in the regions appropriate to each Q 2 bin by the form x~ x . Due to the simplified form of the 
HI parametrization and the extensive x range of the data fitted (which drifts to larger x as Q 2 
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increases), strictly speaking these values should not be compared with Ag. Nevertheless they 
do give a good representation of the errors as a function of Q 2 . 

We saw in Table |] that the x 2 values for the R 2j 4 descriptions of the BCDMS data 
are significantly worse than those for Ri^. The reason, as already mentioned, is that it is 
these data that constrain a s (M|) to be 0.113 in the global fit. The sets with a s (M|) = 0.120 
lead to scaling violations for the structure function which are stronger than those measured 
by BCDMS. This is illustrated in Fig. 4, which compares the Ri and R2 descriptions of the 
"medium 2" SLAC [^T| and BCDMS [|Kj proton structure function measurements. Only those 



data points which pass the W 2 > 10 GeV 2 cut (and which are therefore included in the global 
fit) are shown. The value of a s is mainly constrained by the BCDMS data (on both hydrogen 
and deuterium targets) in the x = 0.35, 0.45 bins. Here the preference of the data for the Ri 
set is evident. Although the SLAC data have a slight preference overall for the larger a s value 
(see Table H), the error bars are larger and so the statistical significance is weaker. There is also 
the possibility that the lower Q 2 data are contaminated with higher-twist contributions, which 
would of course distort the a s measurement. However, a careful analysis || of the combined 
SLAC and BCDMS data, incorporating a phenomenological higher-twist contribution, yields 
the value 

a s {M%) = 0.113 ±0.005 , (8) 

where the statistical and systematic errors have been combined. 

The NLO QCD fits to the medium/large x fixed-target structure function data are relatively 
insensitive to the gluon distribution which, in the MRS global fits, is constrained in this x 
region by the large pt prompt photon data, see for example Ref. p4j| . However at small x 



there is a significant a s - gluon correlation, since dF2/d In Q 2 ~ a s P q9 <8> g- This is evident in 
Fig. 5(a), which shows the gluon distributions of the four R fits in the small x HERA region at 
Q 2 = 5 GeV 2 . At small x the Ri gluon is evidently some 10 — 20% larger than the R2 gluon. 
The same hierarchy is visible in the R3 and R4 gluons, but these are both slightly steeper. 
The reason is that constraining X g to be equal to Ag forces a larger value for the former, see 
Fig. 3, and a correspondingly steeper gluon. Given that the HERA data clearly favour different 
values for X g and Ag (Table we may take the spread in the Ri, R 2 gluons to represent the 
uncertainty in the distribution at small x. It is interesting to note that the previous MRS (A') 
gluon lies within this band. 

Also shown in Fig. 5(a) is the GRV "dynamical" gluon ||. This is much larger than any 
of the four R gluons, and therefore yields a more rapid Q 2 evolution for the structure function. 
This was already evident in Fig. 2, which also showed that such a strong Q 2 dependence is now 
ruled out by the HI and ZEUS data. 

Figure 5(b) shows the relative behaviour of the four R gluons at Q 2 = 10 4 GeV 2 , the 
scale relevant for the Fermilab jet data. As expected, the curves have evolved closer together, 
although the qualitative differences seen at Q 2 = 5 GeV 2 are still apparent. Quantitatively, we 
see that the crossover points of the Ri and R3 gluons (with the same a s ) have moved to smaller 
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x in accordance with the evolution equation, dgjdXwQ 2 ~ P gg (g> g, where g is sampled in the 
convolution at slightly larger x. 

Figure 6 compares the NLO QCD predictions^] for the single jet inclusive Et distribution 
at the Tevatron, using the new Ri and R2 sets, with data from the CDF and DO || col- 
laborations. Only the statistical errors on the data are shown, and in each case the data have 
been renormalized by a constant factor to give the best overall fit.Q As the perturbative QCD 
predictions become unstable at small Et, only data with Et > 50 GeV are included in the 
comparison. The corresponding x 2 values are listed in Table 0. The first and most obvious 
point to note is that highest E T CDF points tend to lie above the theoretical prediction, a 
tendency which is not evident in the DO data. The "disagreement" with the CDF data has 
been the subject of much recent discussion, explanations ranging over a variety of new physics 
effects |26| (but see also [|14], ^7\, p5|). In view of the fact that the large E T DO data appear 
perfectly compatible with the R12 predictions, we believe that it is premature to draw any firm 
conclusions about disagreements with perturbative QCD. More interesting is the discriminating 
power of the high-precision Et < 200 GeV jet data. Since the Q 2 evolution length is so long, the 
steepness of the parton distributions in x, and therefore of the jet Et distribution, is strongly 
correlated with the value of a s |T| . The effect is clearly seen in the differences between the Ri 
(a s (M|) = 0.113) and R 2 (a s (M|) = 0.120) predictions in Fig. 6. One could therefore, at least 
in principle, use these data to measure a s , see for example Ref. [|13|]. Such an analysis would 
however require a careful consideration of the systematic errors on the shape of the measured 
Et distribution, and also of the correlation with the gluon distribution, which is responsible 
for a significant fraction of the uncertainty in the cross section in this Et range. For instance, 
if we consider the Ri and R3 descriptions (which have the same a s ) of the medium- Et CDF jet 
data, the reduction in \ 2 from 222 to 169 in Table |2| is due primarily to the difference in slope 
of the gluons in the x ~ 0.1 — 0.2 region, as shown in Fig. 5(b). However the slope of the jet 
cross section is more sensitive to a s than to the change in the gluon distribution, as reflected 
in the reduction of x 2 to 52 for these jet data for parton set R2. 

Taken at face value, the CDF and DO jet data would appear to favour a larger value of a s , 
of order 0.122 and 0.118 respectively. This suggests that a global analysis incorporating these 
data would yield a set of partons very similar to R 2 with a s (M 2 z ) « 0.120. 

In conclusion, we are able to obtain an excellent and economical NLO DGLAP-based de- 
scription of a wide range of deep inelastic and hard scattering data over an increasingly large 
kinematic domain, namely down to Q 2 = 1.5 GeV 2 and x ~ 10~ 4 . The new HERA data have 
considerably improved our knowledge of the gluon in the small x region. We present^ two sets of 
partons Ri and R 2 corresponding to two values of a s , respectively a s (M|) = 0.113 (the "classic" 

4 The calculation uses the next-to-leading-order parton level Monte Carlo JETRAD and the cuts and jet 
algorithm applied directly to the partons are modelled as closely as possible to the experimental set-up. We 
thank Nigel Glover for help in performing these calculations. 

5 Note that the normalization factors shown in Fig. 6 are in each case well within the quoted experimental 
uncertainty. 

6 The FORTRAN code for the R sets is available by electronic mail from W.J.Stirling@durham.ac.uk 
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DIS value determined by the scaling violations of the fixed-target data), and a s (M|) = 0.120 
(favoured by the HERA F 2 and Fermilab jet data). It is clear that including the jet data, with 
statistical errors only, in the global analysis would in fact discriminate between our two sets[] 
of partons in favour of R2. 
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Figure Captions 



Fig. 1 The description of the EMC and preliminary HI |21J data for F£(x, Q 2 ) by the Ri 
set of partons. The other Rj sets of partons give essentially identical values of F£. 

Fig. 2 The continuous and dashed curves correspond to the values of the proton structure func- 
tion F 2 obtained from the Ri and R2 sets of partons (which have, respectively, QCD 
couplings corresponding to ot a (M|) = 0.113 and 0.120) at twelve values of x chosen to be 
the most appropriate for the new HERA data. For display purposes we add 0.5(12 — i) to 
F2 each time the value of x is decreased, where i = 1, 12. For comparison the dotted curves 
show the prediction obtained from the GRV set of partons |J. The experimental data 
are assigned to the x value which is closest to the experimental x bin. Thus the ZEUS 
data[|J are shown in groupings with x values (3.5,6.3,6.5 x 10 -5 ), (1.02,1.20 x 10~ 4 ), 
(1.98, 2.53 x 10- 4 ), (4.0, 4.5 x 10~ 4 ), (6.32, 8.00 x KT 4 ), (1.02, 1.20 x 10" 3 ), (1.612 x lO" 3 ), 
(2.53,2.60 x 10~ 3 ), (4.00 x lO" 3 ), (6.325 x 10~ 3 ), (1.02 x 10" 2 ), (1.612 x 10~ 2 ), and the 
HI data § in groupings with x values (3.2,5.0 x 10~ 5 ), (0.80,1.30 x 10" 4 ), (2.0,2.5 x 
10" 4 ), (5.0 x 10- 4 ), (6.3,8.0 x KT 4 ), (1.3 x 10~ 3 ), (1.585 x 10~ 3 ), (2.0,2.5,3.2 x 10~ 3 ), 
(3.98,4.0,5.0 x KT 3 ), (6.3 x KT 3 ), (8.0 x KT 3 ), (1.3 x KT 2 ). The E665 data f22[, which 
are shown on the curves with the five largest x values, are measured at x — (2.46 x 10~ 3 ), 
(3.698, 5.2 x 10- 3 ), (6.934 x 10" 3 ), (8.933 x 10~ 3 ), (1.225, 1.73 x KT 2 ). 

Fig. 3 The exponents X g and Xs as a function of Q 2 calculated from the R2 (upper plot) and 
R4 (lower plot) sets of partons. Also shown are the values of A obtained by the HI 
collaboration || by fitting their data to the form F2 = Ax~ x for x < 0.1 at different 
values of Q 2 . 



Fig. 4 The fit to the "medium" x SLAC Q and BCDMS data for Ff by the R x and R 2 sets 
of partons which have, respectively, QCD couplings corresponding to a s (M§) = 0.113 
and 0.120. 

Fig. 5 The gluon distribution at (a) Q 2 = 5 GeV 2 (upper figure) and (b) 10 4 GeV 2 (lower figure). 
The bold curves correspond to the new R 1( R 2 , R3 and R 4 sets of partons. For comparison 
we also show the MRS(A') E21 and GRV M gluons for x < 10" 3 at Q 2 = 5 GeV 2 . 



Fig. 6 The next-to-leading-order QCD description of the CDF J7| and DO single jet inclusive 
Et distribution by the Ri and R 2 sets of partons. These data are not included in the 
global analysis. The overall normalization of the QCD predictions is fitted to the data 
sets and the value found in each case is shown on the individual plots. 



10 




Fig. 2 



0.3 



F p 



0.2 



1 1 

1 


i i i i i 
x=0.35 










MRS(Rl) 
MRS(R2) 

i 








10 



Q 2 (GeV 2 ) 



10 



0.2 



F p 



0.1 



0.1 



F p 



0. 05 



1 III 
1 III 


i i i i i i i i i 
x=0.45 

* • • 

I I I I I I I I 


10 Q 2 (GeV 2 ) 1q2 


i iii 



1 1 1 1 1 1 1 1 1 

x=0.55 
• •• * • • • . # # 

i i i i 



10 



Q 2 (GeV 2 ) 



10 



0.04 



0.02 



1 1 1 1 H^l >. N 


i i i i i i i i 
^ x=0.65 


i 








£ SLAC 






+ BCDMS(x0.98) 




i i 


i 


10 


2 

Q 2 (GeV 2 ) 10 


Fig. 4 







20 

xg(x) 

15 



(a) 



10 



MRS(A') 



MRS(Rl) 
MRS(R2) 
MRS(R3) 
MRS(R4) 



Q 2 =5 GeV 2 



10 



.-4 



1.1 



(b) 



10 3 X 10 2 



10 



Q 2 =10 4 GeV 2 



-l 



1.05 



s 

^ 1 

s 

0. 95 



Rl 



R3 



R2 



/ 

/ 

>/ 



0.9 



10 



.-4 



R4 



10 



-3 



10 



-2 



10 



-1 



Fig. 5 



1 .4 i 

1 .2 - 

MRS(Rl) 
1 _ and CDF data (x0.914) 

0.8 - 
0.6 - 



0.4 
0.2 


0.2 



1 .4 i 

1 .2 - 

MRS(R2) 
1 _ and CDF data (x 1.008) 

0.8 - 



0.6 

0.4 







0.2 - 



-0.2 - 



-0.4 



1 .4 



100 200 300 400 
E T (GeV) 



-0.4 



1 .4 



100 200 300 400 
E T (GeV) 



.2 - 

MRS(Rl) 
1 _ and DO data (x0.906) 



1 .2 - 

MRS(R2) 
1 _ and DO data (x 1.000) 



0.8 
0.6 
0.4 
0.2 

- 

-0.2 -- 
-0.4 



100 200 300 400 
E T (GeV) 



0.8 - 
0.6 
0.4 
0.2 

- 
-0.2 - 

-0.4 



100 200 300 400 

E T (GeV) Yig. 



